Dynamics of charged gibbsite platelets in the isotropic phase 



Dzina Kleshchanok/ Marco Hcinen,^ Gerhard Nagele,^ and Peter Holmqvist^ 

^ Van 't Ho ff Laboratory for Physical and Colloid Chemistry Debye Research Institute Utrecht University H.R. Kruytgebouw, 

N-701 Padualaan 8, 3584 CH Utrecht, The Netherlands 
^Institute of Complex Systems (ICS-3), Forschungszentrum Jillich, D-52425 Jiilich, Germany 

(Dated: September 16, 2011) 

We report on depolarized and non-depolarized dynamic light scattering, static light scattering, 
and static viscosity measurements on interacting charged gibbsite platelets suspended in dimethyl 
sulfoxide (DMSO). The average collective and (long-time) translational self-diffusion coefBcients, and 
the rotational diffusion coefRcient, have been measured as functions of the platelet volume fraction 
(j), up to the isotropic-liquid crystal (I/LC) transition. The non-depolarized intensity autocorrela- 
tion function, measured at low scattering wavenumbers, consists of a fast and a slowly decaying 
mode which we interpret as the orientationally averaged collective and translational self-diffusion 
coefRcients, respectively. Both the rotational and the long-time self-diffusion coefHcients decrease 
very strongly, by more than two orders of magnitude, in going from the very dilute limit to the I/LC 
transition concentration. A similarly strong decrease, with increasing (p, is observed for the inverse 
zero-strain limiting static shear viscosity. With increasing increasingly strong shear-thinning is 
observed, accompanied by a shrinking of the low shear-rate Newtonian plateau. The measured diffu- 
sion coefRcients are interpreted theoretically in terms of a simple model of effective charged spheres 
interacting by a screened Coulomb potential, with hydrodynamic interactions included. The disk- 
like particle shape, and the measured particle radius and thickness polydispersities, enter into the 
model calculations via the scattering amplitudes. The interaction-induced enhancement of the col- 
lective diffusion coefRcient by more than a factor of 20 at larger (p is well captured in the effective 
sphere model, whereas the strong declines both of the experimental translational and rotational 
self-diffusion coefRcients are underestimated. 



I. INTRODUCTION 

Colloidal platelets are abundant in nature {e.g., as clay 
minerals or red blood cells) and can be readily synthe- 
sized in the laboratory in form of mixed metal hydrox- 
ides, gibbsite, gold platelets, et cetera. The most promi- 
nent examples of colloidal platelets are various types of 
natural clays [1], which figure in sediment transport in 
rivers, and in the oceans and lakes, and which are re- 
sponsible for dangerous landslides [2]. Clays are widely 
used as rheological modifiers for surface coatings, paints, 
and drilling fluids [3, 4]. 

The rheological applications of clays are based on their 
microstructural properties, originating from a highly 
anisotropic shape and a correspondingly anisotropic par- 
ticle interactions. As a necessary step towards under- 
standing the more complex behavior of concentrated clay 
platelet suspensions regarding rheology, sedimentation, 
and sediment transport, the fundamental dynamic prop- 
erties in dilute isotropic solutions must be addressed. 

So far, most studies on the dynamics of clay and col- 
loidal platelet suspensions have been focused on the non- 
equilibrium ergodic to non-ergodic transition in dense 
systems [5-7], and on the properties of the non-ergodic 
state [8, 9]. Only few theoretical and experimental 
investigations have been done so far on less concen- 
trated, fluid-state isotropic suspensions of clays or col- 
loidal platelets. For instance, both the translational and 
rotational diffusion coefficients of a single platelet were 
studied by simulation, and the resulting data for an ex- 
tended range of aspect ratios were fitted to general poly- 



nomial expressions for cylinders by Ortega and Garcia de 
la Torre [10]. This study allows to compare the single- 
platelet diffusion properties of non-interacting platelets 
with the results of experimental diffusion measurements, 
but it includes no hint on how diffusion is affected by 
particle interactions. 

The concentration dependence of sedimentation and 
diffusion coefficients of uncharged, sterically interacting 
platelets and clays has been explored in [11, 12] by 
dynamic light scattering and ultracentrifugation. For 
platelets, the earlier investigations in [11] were made in 
a rather dilute concentration regime where no effect of 
particle interactions on diffusion was detected. More- 
over, explicit theoretical expressions have been derived 
for the time-dependent intensity autocorrelation function 
of non-interacting cylindrical particles such as platelets 
[11, 13-16] which, in principle, can be used to deter- 
mine the single-platelet rotational diffusion coefficient 
from standard dynamic light scattering or X-ray photon 
correlation spectroscopy (XPCS) measurements in non- 
interacting particle systems. 

For rather dilute systems of uncharged polymer- 
grafted clay suspensions [12], normalized intensity au- 
tocorrelation functions (lACFs) have been measured in 
vertical-vertical (VV) and vertical-horizontal (VH) scat- 
tering geometry [17]. A strong slowing of the diffusive 
modes was found both in VH and in VV geometry. In 
VV geometry, only a single diffusive mode was detected, 
even though it was argued that a second, cooperative 
mode should be present due to the osmotic pressure of 
polymers grafted on the clay particles [12]. 

In the present study, we explore how the dynam- 
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ics of charged colloidal gibbsitc platelets suspended in 
DMSO and present in their isotropic phase, is affected 
by direct and hydrodynamic particle interactions. Using 
static light scattering (SLS) and depolarized and non- 
depolarized dynamic light scattering (DDLS and DLS), 
we study the effect of the platelet concentration on the 
measured translational and rotational self-diffusion coef- 
ficients, and on the collective diffusion coefficient. More- 
over, the concentration- and shear-rate dependence of the 
static dispersion viscosity is determined experimentally. 
Our scattering and rheological experiments cover the full 
isotropic phase concentration regime up to the I/LC tran- 
sition. 

The theoretical description of the platelet dynamics 
in non-dilute systems is severely complicated by the 
anisotropy in the direct [IS] and, to an even larger ex- 
tent, by the indirect hydrodynamic interactions (His). 
No readily applicable theoretical schemes or simulation 
results are available for our gibbsite system. There- 
fore, in a first attempt to analyze theoretically the ex- 
perimental trends in the diffusion data, we use a sim- 
plifying model where the direct interactions in a pair 
of gibbsite platelets are described approximately by the 
spherically symmetric, repulsive part of the Derjaguin- 
Landau-Verwey-Overbeek (DLVO) potential [19]. Clays, 
such as Laponite, are mixed silica and metal oxides. 
This can lead to charges of opposite sign on their faces 
and edges [20-22]. The present gibbsite platelets con- 
sist solely of aluminium hydroxide, so that their faces 
and rims of gibbsite have surface charges of the same 
sign [23]. In our effective sphere model of gibbsite in 
DMSO, we account for the influence of the His but leave 
out the effects of interaction polydispersity. However, 
the disk-like shape, and the experimentally determined 
size polydispersity of the gibbsite platelets, are accounted 
for in the scattering amplitudes. The simplifying model 
of cylindrical platelets with spherically symmetric direct 
and hydrodynamic interactions permits us to take ad- 
vantage of a wealth of theoretical and computer sim- 
ulation results on rotational [24, 2-'')] and translational 
self-diffusion, and collective diffusion of charged colloidal 
spheres [25-27]. The model can be expected to apply 
approximately to semi-dilute suspensions of isotropically 
dispersed platelets well below the overlap volume frac- 
tion. For larger concentrations, it is bound to fail. 



II. EXPERIMENTAL DETAILS 

We study suspensions of hexagonal colloidal gibbsite 
[7-Al(OH)3] platelets of average radius (R) = 44.2 nm 
and average thickness (h) ~ 7.66 nm, dispersed in 
dimethyl sulfoxide (DMSO). The gibbsite particles can be 
described approximately as a polydisperse system of thin 
circular cylinders of radius R and (small) height h. We 
have determined the radius and height distributions from 
electron microscopy pictures, and using atomic force mi- 
croscopy, giving polydispersities (i.e., relative standard 



deviations) in R and h of Sfi = 17.3% and Sh = 55.3%, 
respectively. 

The mean aspect ratio, p = (h) /(2 (R)) = 0.087, 
of platelets is so small that the ultrathin disk limit 
(h — > 0) can be applied as a reasonable approximation. 

In this hmit, one obtains = (D°'" + 2D°-^)/3 = 
kBT/{12r]o {R}) « 3.8 x lO^^^ ^^2/3 for the orientation- 
ally averaged single-disk translational diffusion coefficient 
[28, 29], and D^'-^ = ikeT / {i2r^o{Rf) w 2.2 x lO^/s for 
the (end-over-end tumbling) rotational diffusion coeffi- 
cient determined in small-g DDLS experiments on non- 
interacting platelets. Here, ks is the Boltzmann con- 
stant, T is the absolute temperature, and 770 is the solvent 
shear-viscosity. For p — )■ 0, D^'^ becomes equal to the 

rotational diffusion coefficient, -Dr", which characterizes 
rotation with respect to the platelet rotational symme- 
try axis. A small but finite value of p lowers somewhat 
the values of D° and -D°'^, and D^'" and D"'^ become 
different from each other. Tirado and Garcia de la Torre 
[30, 31] provide precise polynomial fits to their simulation 
data for the diffusion coefficients of cylindrical platelets 
as a function of p. However, the aspect ratio range, 
p > 0.1, covered by most of these fits, does not include 
the small aspect ratio of the present gibbsite platelets. 
We can alternatively estimate the effect of a non-zero 
p using analytic results [32] for the single-particle diffu- 
sion coefficients of oblate spheroids of same p and same 
volume Vp = Ti^R^h), giving D^'-"- = 1.4 x 10^/s and 
DO = 3.1 X 10-i2mVs. 

DMSO is a polar, aprotic, and low-viscous solvent of 
dielectric constant e = 47.2, at T = 293i^, and viscosity 
770 = 2 X 10^"^ Pa-s, in which the platelets are charge- 
stabilized, forming a suspension that remains transpar- 
ent up to the isotropic-liquid crystal (I/LC) transition. 
Different from aqueous gibbsite suspensions, which are 
turbid already at low concentrations, gibbsite in DMSO 
systems are transparent and show no multiple scattering 
in the investigated concentration range. Moreover, differ- 
ent from aqueous suspensions, there is no residual CO2 
contamination, and no solvent self-dissociation, so that 
low-ionic strength systems can be easily prepared. Thus, 
we can use DDLS to study the rotational and transla- 
tional diffusion as a function of the gibbsite number con- 
centration up to the I/LC transition point. Our standard 
(D)DLS/SLS apparatus is equipped with a krypton ion 
laser of wavelength Aq ~ 647 nm as a light source, and a 
A/2 plate used as polarizer and analyzer (Bernhard Halle 
Nachfl., Berlin, Germany). Each sample was measured 
both in vertical-unpolarized (VU) and VH scattering ge- 
ometry, for values of the wave number, g, smaller than 
the value, q^, where the primary peak of the mean scat- 
tered intensity, /(?), occurs. 

DLS and static viscosity data were recorded for a large 
number of gibbsite volume fractions 4> = Af/(^?m^) up to 
the I/LC transition at 0i/lc ~ 8%. Here, M is the total 
mass of added gibbsite of known mass density g^m and 
V is the suspension volume. In addition, the static shear 
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viscosity, r], was measured as a function of shear rate 7, 
for a large number of concentrations. 



III. RESULTS 
A. Dynamic light scattering results 

The dynamic hght scattering functions in VH and VU 
geometry have been recorded for a large number of scat- 
tering wavenumbers and concentrations up to the vol- 
ume fraction at the I/LC transition. Results for the 
normalized intensity autocorrelation functions (lACFs), 
g2^{q,t), (filled symbols) and g2^{q,t) (open symbols), 
for two different volume fractions = 0.8% and 4.8%, 
are shown in Fig. 1. We point out here that w 
to excellent accuracy, due to the very small optical 
anisotropy of the gibbsitc platelets. The rather dilute 
system at cj) = 0.8% has only weak inter-platelet corre- 
lations, whereas the system a.t <f> = 4.8% is not very far 
from the I/LC transition. 

Consider first the low-concentrated system in panel (a) 
of Fig. 1, where (g, t)~l and gY^{q, t) — l are shown 
for four different 5- values located well below with 
q (R) = 1.26 for the largest considered q. At this low 
concentration, both lACFs are only slightly stretched, 
decaying basically single-exponentially. The stretching is 
due here to the rather small polydispersity, sr = 17.3%, 
in the disk radius R. Due to the small value of p, the more 
pronounced polydispersity in h of value Sh = 55.3% has 
only little influence on the particle diffusion coefficients, 
and on the scattering amplitudes in the experimentally 
accessed g-range where q (R) < 1 and thus q (h) < 0.1. 
The time relaxation of g2^{q,t) is faster than that of 

To globally account for polydispersity and particle cor- 
relation effects, which both give rise to a potentially con- 
tinuous superposition of correlation times, we have fit- 
ted all of the measured functions 32(9, ^) — 1 in panel 
(a) of Fig. 1 by the Kohlrausch-Wilhams- Watts (KWW) 

stretched exponential form, {exp[— (t/ri) ]} , character- 
ized by the overall decay times Ti and stretching ex- 
ponents f3i < 1. An lACF decaying nearly single- 
exponentially has a value of /3 close to one. For example, 
one finds /3i « 0.9 for the lACF's of the low-concentration 
systems in panel (a) of Fig. 1. The KWW mean relax- 
ation rates F^, with i S {VH,VU}, are obtained from 
Ti ~ 1/ [Ti X r{l/ f3i) / f3i\, where F denotes the gamma 
function[.').:)]. In all our experiments, the area coherence 
factor C in the Siegert relation Cgl{q,t) = g2{q,t) — 1 is 
practically equal to one. 

In panel (b) of Fig. 1, the resulting relaxation rates 
FvH and Fyu at </) = 0.8% are plotted as functions of q"^. 
At the small g- values considered where q{R) < 1, the 
relaxation rates in VH and VU«VV geometry can be ex- 
pressed as FvH = QD;^ + O(g^) and Fyy = A?^ + ©(g^), 
respectively [13, 15-17, 34], where Dt and Z?^ are inter- 
preted as the translational and (end-over-end tumbling) 



rotational diffusion coefficients, respectively. In this in- 
terpretation, polydispersity effects are accounted for in 
an overall way through the KWW relaxation rates. Con- 
tributions to Fvv of 0{q'^) which arise from rotational 
diffusion and rotational-translational coupling, are of no 
relevance for the considered small (/-values. 

Regarding the low-concentration system results in pan- 
els (a) and (b), no distinction is required between short- 
and long-time diffusion coefficients. The slopes and in- 
tercepts of Fvu and Fvh in panel (b), obtained from a 
q^ ^ extrapolation, are approximately equal to the 
zero-concentration diffusion coefficients. A zeio-cj) ex- 
trapolation of all investigated systems, with the relax- 
ation rates obtained as described above, leads to D'} = 
2.1 X 10-12 mVs and ZJO-^ = 0.34 x lOVs. Both val- 
ues are somewhat smaller than those obtained from the 
ultrathin platelet approximation using R ~ (R). 

The data for Fyn in panel (b) show an overall linear 
increase in g^, but are more noisy than those for Fyu- 
This originates from the much lower scattering intensity 
in VH geometry, which for the gibbsite system is by a 
factor of 10"^ smaller than that in VU geometry. 

At larger platelet concentrations, where direct and 
hydrodynamic interactions come into play, one notices 
interesting changes in the scattering functions. The 
normalized VU and VH lACFs, and the corresponding 
KWW relaxation rates at = 4.8%, are depicted in pan- 
els (c) and (d) of Fig. 1, respectively. It is apparent from 
panel (c) that, in VU geometry, an additional quickly 
relaxing mode occurs. Moreover, the decay of the slow- 
mode in g^^iq,t) — 1 is slightly more stretched. In VH 
geometry, no additional mode is seen at this larger con- 
centration, but the decay is more stretched than in the 
(f) = 0.8% case. These features arc indicators of pro- 
nounced platelet correlations at = 4.8%, requiring now 
the distinction between short- and long-time diffusion 
properties. The subdiffusive transition regime is char- 
acterized by the structural relaxation time defined by 



1/(6D0 



0.05 ms. 



Using again the Siegert relation, we can fit the two- 
mode decay of observed in panel (c) in the probed 
time window by the two-exponential form. 



1 = 



Ae 



(1) 



involving three fit parameters A, F^^ and F^^ with 
< ^ < 1 and Fyu > F/" > 0. Both the fast-mode 
and the slow-mode relaxation rates, F^^ = q^Dc and 
TJ^ = q^Dt, show the expected diffusive g^-depcndence, 
allowing for the determination of the associated diffusion 
coefficients Dc and Dt, respectively. For (f> = 4.8%, the 
ordering relations Dc > D^ and Dt < D^ are obeyed. In 
panel (d), F^" - = q^jj^ + 0(^4^ plotted as a 
function of q^, with Dt inferred from the small-g slope. 
The corresponding fast mode rate, F^^, is not shown in 
the figure. Its associated diffusion coefficient, D^ has 
been determined as the slope of r^^{q'^), extrapolated 
to g = 0. 
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Figure 1. Panels a and c: Reduced intensity autocorrelation function, g2{q,t) — 1, versus correlation time t, for gibbsite in 
DMSO. Panels b and d: (mean) relaxation rates, Fvh and Tvu, as functions of q^. The gibbsite volume fraction is = 0.8% 
in panels a and b, and (j> = 4.8% in panels c and d. The wavenumbers in a and c are: q = 7.6 x 10~^ nm~^ (black circles), 
q = 1.47 X 10"'^ nm~^ (red squares), q — 2.08 x 10~^ nm~^ (green triangles), and q = 2.84 x 10~^ nm~^ (blue triangles). Open 
symbols: VH scattering geometry results. Filled symbols: VU ~ VV scattering geometry results. 



Panel (d) shows additionally the relaxation rate, 
= 6D:jr+Dtq^ , of the somewhat stretched out single- 
exponential decay of gj^{q,t), determined again using 
the KWW analysis. Note here that the small-(j slope of 
r^^(g^) is equal, within the experimental noise, to that 
of rY^{q^). The slope in both geometries is identified as 
the translational self-diffusion coefficient, Dt, of interact- 
ing platelets in the isotropic phase. Likewise, the coeffi- 
cient D:^ determined from the zero-g intercept ofr^"(g2) 
in panel (d), obeying D:jr < D^'-^ = D^[4> 0), can be 
interpreted as the rotational self-diffusion coefficient at 
non-zero concentrations. 

The same evaluation procedure as explicated in Fig. 1 
for two specific concentrations, was applied to an ex- 
tended set of concentrations up to the I/LC transi- 
tion value. The obtained reduced diffusion coefficients, 
A/-D?> D^/D°-^ and DJD'^, are depicted in Fig. 2 in 
their dependence on 0. According to Fig. 2, Dt and 
D:^ remain constant, within the experimental scatter, 
in the lower-concentration range of < 0.5%. This 
is the expected behavior of particles at low concentra- 
tion which diffuse practically uncorrelated. On the other 
hand, already at (/) = 1%, the fast-mode coefficient, Dc, 
exceeds largely its zero-concentration value Z?^ whereas, 
as viewed on the extended vertical scale of Fig. 2, Dt and 
D:^ are just about to start their strong decline below their 
respective infinite dilution values. 

A concentration dependence similar to that of Dc has 
been reported for the collective diffusion coefficient of 
polymers and flexible rod systems [3">, 3()]. However, 



the gibbsite platelets studied in the present work include 
no flexible parts. Moreover, the onset of the strong in- 
crease of Dc{4)) is observed at concentrations far below 
the platelet overlap concentration, 0* = (3/2)p, of about 
13%. Incidentally, an alternative definition of the over- 
lap concentration invoking the random close packing vol- 
ume fraction, / = 0.64, of monodisperse spheres given by 
£»*(47r(i?)V3) = / results in 0* = (3 {h) /4 {R))f w 8.3%. 
This value is somewhat fortuitously close to the concen- 
tration, 0I/LC = 8%, where the I/LC transition is ob- 
served experimentally. 

Since the dynamics of freely rotating charged platelets 
at concentrations well below the overlap concentration 
(j)* is most strongly influenced by the monopolar terms 
in the far-field electrostatic and hydrodynamic interac- 
tions, and since collective diffusion is to a large extent 
triggered by the osmotic compressibility, Dc{(f>) can be 
expected to behave similar to the collective diffusion co- 
efficient measured in low-salinity suspensions of spherical 
charged colloids [2G, 27]. The latter coefficient also shows 
a distinct rise at lower 0. 

We argue that Fig. 2 depicts the long-time values of 
the translational self-diffusion coefficient Dt which, for 
interacting systems, can be substantially smaller than its 
short-time counterpart, whose deviations from the infi- 
nite dilution limit D^ are comparatively small. That Dt 
should be interpreted as a long-time property indeed fol- 
lows from Eq. (1) and the considered small-g range in 
combination with the observation that Dc > D^ > Dt 
s-t ^ 0.5%. Due to the fast decay of the collective 
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Figure 2. Volume fraction dependence of various dynamic 
properties including the normalized fast-mode diffusion coef- 
ficient, DcjDl, (filled squares), and the slow-mode diffusion 
coefficient, Dt/Dt, (filled circles), both obtained in VU ge- 
ometry, and the normalized (mean) rotational diffusion co- 
efficient, /D^'^ (open squares), obtained from VH mea- 
surements. Stars: Inverse low-shear-rate viscosity (rj/'qo)"^, 
obtained from rheometric measurements. Solid black and 
dashed red lines: Predictions for Dc/Dt and Df /Dt within 
the effective charged sphere model, calculated using the cor- 
rected ($7 scheme. Solid blue line (inset only): Scaling relation 
Df /D^ = I — 1.3(/>^ff for the short-time rotational diffusion co- 
efficient of effective charged spheres, using (^eff = 1.930. Solid 
red line: Df /Df of effective charged spheres, calculated using 
the simplified MCT scheme with far-field His correction. The 
inset magnifies the details at low on a double linear scale. 



mode, is determined essentially by the coefficient 

Dt for times t > l/{q^Dc), which is part of the long-time 
regime t 3> ~ 0.05 ms once the zcro-g extrapolation 
is made in determining Df. Different from self-diffusion, 
the long-time value of is practically equal to the short- 
time value, for all concentrations in the isotropic liquid 
state regime [37, 38]. 



B. Simplifying difFusion model of platelets 

To analyze the trends in the concentration dependence 
of the gibbsite platelet diffusion coefficients, and to gain 
further support for the interpretation of Dc as a collec- 
tive difFusion coefficient, we introduce here a simplifying 
model. In this model, the platelets are described regard- 
ing their interactions as uniformly charged effective hard 
spheres of radius a, interacting by the repulsive electro- 
static part, 




/3Mei(r- > 2a) = LbZ"^ 



1 + Ka 



(2) 



q / nm 



Figure 3. Normalized mean scattered intensity, I{q) / [locj)), 
of gibbsite platelets in DMSO, for (f) = 0.16% (black squares) 
and (j) = 0.88% (red circles), in comparison to the theoretical 
fits using (/)eff = 0.43% (black line) and <j)cs — 1.68% (red line). 
Black dotted curve: measurable form factor Pm{q)- Blue 
dashed curve: decoupling amplitude X{q). Insets: experi- 
mental histograms of the platelet radius R (top) and thick- 
ness h (bottom). The red shaded areas are least-squares fits 
by unimodal Schulz distributions with mean values (R) = 44.2 
nm and (h) = 7.66 nm, and relative standard deviations of 
sr — 17.3% and Sh ~ 55.3%, respectively. 



of the DLVO potential [I:/]. Here Lb = 1-18 nm is the 
Bjerrum length of DMSO at T = 293X, Z is the effective 
particle charge in units of the proton elementary charge, 
and = 4:TrLB[ncs\Z\ + 2ns] /{I — 4'cs) is the square of 
the Debye screening parameter. The latter is determined 
by the number concentration, n^, of residual 1-1 elec- 
trolyte ions which we assume to be equal to and the 
number concentration, ncs\Z\, of surface-released mono- 
valent counterions. Furthermore, Ues = 30cff/(47ra^) is 
the number concentration of effective spheres. There are 
several alternatives to define the effective sphere radius 
a, depending on the quantity considered. In the present 
model, a is obtained from equating the 2'^'* virial coeffi- 
cient of neutral spheres of radius a to that of cylindrical 
platelets of radius {R) and height (h), assuming all ori- 
entations to be equally probable [30]. This results in 
a = 0.735 (R) = 32.5 nm. A slightly smaller value of 

1/3 

a = (37r/32) ' (R) ~ 29.4 nm would be obtained from 
the 2"^^ virial coefficient in the zeio-h limit. 

In addition to the DLS measurements described before, 
we have measured the static mean scattered intensity, 
/(g), using SLS. Fig. 3 includes the experimentally de- 
termined intensities (open symbols), for platelet volume 
fractions cj) = 0.16% (black) and 0.88% (red). 

Light scattering measurements are restricted to a 
smaller g-range than small angle X-ray scattering (SAXS) 
so that the interaction peak in I(q) is resolved experi- 
mentally for the lower concentrated system only. Using 
the effective sphere interaction model, with Uei{r) accord- 
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ing to Eq. (2), we have calculated I{q) approximately 
within the rotational-translational decoupling approxi- 
mation [37, 40], where it is given by 

/(g) -/o0P™(g)5,„(g), (3) 

with the so-called measurable static structure factor, 

Sm{q)^[l^X{q)]+X{q)S{q). (4) 

Here, /q is a ^-independent factor with the dimension of 
intensity, which is the same in all intensity measurements 
corrected for source intensity, recording time, and scat- 
tering volume. Moreover, S{q) is the structure factor of 
the monodisperse system of effective spheres interacting 
by the pair potential in Eq. (2). The measurable form 
factor, Pm{q), and the dimensionless decoupling ampli- 
tude X{q), 

Pmiq) = {f)iq), (5) 
X{q)^{ff{q)/P^{q), (6) 

with < X{q) < 1, are given in terms of the first and sec- 
ond moments of the scattering amplitude, f(q, R, h, ji) = 
Vpb{q, R,h, fi), of a cylindrical platelet of volume Vp = 
TrR^h, where iJ. denotes the cosine of the angle between 
a platelet's rotational symmetry axis and the scattering 
wave vector. The simple analytic expression for the di- 
mensionless form amplitude, b{q,R, h,fi)^ is given in the 
literature [15, 41, 42]. 

Eqs. (3)- (6) are based on the assumption that the par- 
ticle shape, orientation, and the size polydispersity affect 
the form amplitudes only, which excludes anisotropic in- 
teraction effects. Therefore, the decoupling approxima- 
tion can be reasonably applied at lower concentrations 
only, where the platelets rotate essentially freely. The 
brackets, <...>, in Eqs. (5) and (6) represent averages 
with respect to the platelet orientation variable /x (as- 
suming free rotation), and the platelet radius and height 
joint distribution function p(i?, /i), which we assume to be 
factorisable as p{R, h) = Pr{R) x ph{h). Histograms for 
the marginal distributions pi^{R) and Ph{h) of the dried 
specimen, obtained from electron microscopy and atomic 
force microscopy, respectively, are shown in the insets of 
Fig. 3. On fitting unimodal Schulz distributions both to 
Pr{R) and ph{h), we obtain mean values (i?) = 44.2 nm 
and (h) ~ 7.66 nm, at relative standard deviations of 
55.3% and 17.3%, respectively. Note that for q<2/ (R), 
X{q) K, 0.72 stays practically constant. 

The static structure factor S{q) in Eq. (4), describing 
pair correlations of effective charged spheres, has been 
calculated using the modified penetrating-background 
corrected rescalcd mean spherical approximation (MPB- 
RMSA), introduced recently by part of the present au- 
thors [43, 44]. This analytic method allows for an effi- 
cient and accurate calculation of structure factors of non- 
overlapping spheres, interacting by a repulsive Yukawa- 
type potential such as the one given in Eq. (2). 

A fit of the decoupling approximation expression for 
I{q) in Eq. (3) to the experimental intensities involves Z , 



(f>oS, lo and, within a reasonably small range, also Us as 
the fitting parameters. Since an accurate determination 
of the effective particle charge number Z relies on fitting 
the maximum in I{q), not located inside the SLS-resolved 
q-range for most of the considered concentrations, we al- 
ternatively estimate Z using the experimentally observed 
I/LC transition concentration, determined in our low- 
salinity experiments as 0i/lc — 8%. Earlier simulations 
[45, 46] predict the isotropic- nematic transition of neutral 
hard disks to occur at ni/i^QcP « 4, where nj/Lc is the 
isotropic-phase number density at the transition point, 
and d is the disk diameter. Using this relation, we esti- 
mate the charge on a gibbsite platelet by assuming that 
d = 2[{R) +K-1] and 71i/lc = (/)i/Lc/(7r (i?2/i)), with 

= 47ri_B(ni/Lcl-Z'| +2ns)/{l -(/"i/lc)- For the consid- 
ered large concentration at the I/LC transition point, the 
residual salt contribution to k can be neglected. Solving 
for the effective platelet charge gives the value Z = 71. 
For simplicity, and since the concentration dependence of 
Z for gibbsite in DMSO is unknown, this effective charge 
value is used in all our calculations, independent of </>. 
According to [23], Z = 71 is a reasonable charge value 
for gibbsite platelets. 

As demonstrated in Fig. 3, decently good fits of the ex- 
perimental I{q) in the probed g-range are obtained from 
adjusting the only remaining fit parameter, 0cff, entering 
into the MPB-RMSA calculation of S{q). The global fac- 
tor Iq in Eq. (3) only sets the overall intensity scale and is 
therefore system independent. Assuming a homogeneous 
linear relation between 0off and </>, so that 0off = for 
(/) = 0, from our fit we obtain 0cff = 1.93(/). 

While the theoretical fit of the SLS I{q) is satisfying, 
future SAXS measurements are required to obtain I{q) 
in a more extended g-range which covers all its structural 
features, allowing for a more detailed fitting. 

The effective sphere interaction model in combination 
with the rotational-translational decoupling approxima- 
tion is easily generalized from the SLS intensity to the 
normalized time-dependent electric field autocorrelation 
function, which for small q in VV (~VU) geometry is 
given by 

gY^{q,t) = [l-i3]e-9'^W+i?5(q ^ 0)6"'?'^'=^ (7) 

with B = X{0)/[X{0) + (1 - A:(0))5(0)]. Here, W{t) = 
1/6 < [r(i) -r(0)]^ > is the mean-squared displacement 
of a charged effective sphere at center-of-mass position 
r{t) with initial (short-time) slope Df , and final (long- 
time) slope Df, where < Df < D^. Moreover, 
on ignoring the very small difference between long- and 
short-time collective diffusion coefficients, which becomes 
noticeable at very high concentrations only, Dc is given 
by Dc = D^H{q ^ 0)/Siq 0), where H{q) is the 
so-called hydrodynamic function [16, 37, 47]. Eq. (7) is 
fully consistent with the for the gibbsite in DMSO sus- 
pensions observed two- mode decay of the lACFs, fitted 
using Eq. (1), and with the interpretation of the exper- 
imental Dt as a long-time self-diffusion coefficient. We 
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note that in the simplifying model considered here, the 
diffusion coefficients of platelets are simply approximated 
by those of the effective charged spheres. 

Using the effective sphere parameters Ug = 7/iM, a = 
32.5 nm, and Z — 1\ determined as described above, and 
on varying the (effective) volume fraction in small steps 
using the relation 0cff = 1.930, a set of S{qys has been 
generated using the MPB-RMSA scheme, which in turn 
was employed as input to the otherwise parameter-free 
self-part corrected ^7 method of calculating H{q) [2-5, 48- 
5f]. The corrected 5^ method includes many-body hy- 
drodynamic interactions in an approximate way, mak- 
ing predictions for the H{q) of Brownian spheres with 
Yukawa-type repulsion in overall good agreement with 
Stokcsian Dynamics computer simulation results [oU]. 

Regarding the approximations invoked in our simpli- 
fying analytic model, the agreement between the calcu- 
lated Dc{4>) and the experimental data is rather satis- 
fying. The increase of Dc with increasing (jj to values 
about 20 times larger than is qualitatively captured. 
The theoretical Dc reaches a shallow maximum at about 
(j) = 3%, originating from the interplay of osmotic com- 
pressibility and sedimentation coefficient H{q 0) as 
discussed, e.g., in [26, 52]. 

The large-g limiting value of H{q) is equal to the re- 
duced short-time self-diffusion coefficient Df /D^. The 
self-part corrected (57-scheme result for Df{4)) is included 
into Fig. 2 (see also the inset). In the explored (p- 
range, Df{(j)) decreases only mildly with increasing 0. 
Its (/)-dependence is overall well described by the form 
(1 — at(j>^(^)Df with at ~ 2.9, which is the typical 
concentration-dependence of ZJf for low-salinity systems 
of charged spheres [27]. For cj) > 2%, however, the de- 
cay of the experimental (long-time) Dt in Fig. 2 with 
increasing cf) is very strong. For volume fractions (f) > 7% 
near to ^i/lcj Dt has decreased to values less than 1% of 
Dt . This is similar to the slowing down of self-diffusion 
seen for block copolymers in a neutral solvent when the 
disordered-ordered transition is approached [53], in poly- 
mer solutions with increasing concentration [3G], and 
also for polymer grafted clay particles [12]. The mea- 
sured values of Df/Z?° for gibbsite are much smaller 
than those reached by the long-time translational dif- 
fusion coefficient of low-salinity charge-stabilized spheres 
in the fluid phase, which reaches its minimal value of 
Dt /D^ w Dt/Df w 0.1 at the freezing transition vol- 
ume fraction (pf [54, 55]. For the effective sphere model 
used here, 0/ can be estimated on basis of the empirical 
Hansen- Verlet freezing criterion S{qm,4>eS = 4'f) ~ 3.1 
for charged spheres at low salinity [50-58]. With S{q) 
computed in MPB-RMSA, this results in the freezing 
transition volume fraction 0off = 42%, corresponding to 
(j) = 0eff/1.93 = 22%. 

We attribute the strong decay of the experimental Dt 
at larger </) to the uprising influence of the anisotropic 
electro-steric interaction parts, and to the hydrodynamic 
rotational-translational coupling effects in platelet sys- 
tems. These anisotropic interaction effects, not included 



in the effective sphere model, cause an additional strong 
slowing down of the translational and rotational self- 
dynamics. The inset in Fig. 2 displays the experimental 
data for Dt/Dt (and D^ /D^-^) on a smaller concentra- 
tion range cj) < 0.025, in comparison with the corrected 
(57-scheme predictions for Df /D^, and simplified mode- 
coupling theory (MCT) results for Dt /D^ within the ef- 
fective charged sphere model. In the simplified MCT 
result, the enhancing influence of His on Dt, typical of 
low-salinity systems, is accounted for [59-Gl]. Note here 
that the peak height, S{qjn) = 1.7 of the static structure 
factor at (pi/LC = 8% is still rather close to one so that 
our usage of the simplifled MCT solution in place of the 
fully self-consistent MCT solution for Dt is justified. 

The short-time rotational self-diffusion coefficient, D^ , 
in chargcd-sphcrc colloidal systems follows at low salin- 
ity to a good accuracy the scaling relation D^ /D^ = 
1 — ar4>'^ff, with Qr ~ 1.3, in the whole (p-vaiige cov- 
ered in Fig. 2 [24, 25]. This curve is shown as the 
solid blue line in the inset of Fig. 2. Different from 
platelets with /i > 0, spheres are characterized by a 
single zero-concentration rotational diffusion coefficient 
D^ = kBT/{8TTr]oa^). Note here that D^/D^ decreases 
less strongly with increasing (jj than the experimental 
Df/D°. This difference originates from the shorter- 
ranged hydrodynamic self-coupling of rotational motion 
[24]. The experimental diffusion coefficients Dt and D:^ 
depicted in Fig. 2, obey overall the same ordering rela- 
tion D^{<j3) > Dt{4>) as their short-time counterparts in 
the effective sphere model. 

The low-g expression for gi^il^t) in our simplifying 
model is given by (see, e.g., [34, G2]) 

5r('Z,t)=e-'^'^(*)a(i), (8) 

where Grit) =< P2{u.{t) ■ u(0)) > is the rotational self- 
dynamic correlation function of spheres, with the optical 
axis of a sphere characterized by the unit vector u, and 
with P2 denoting the 2"'^-order Legendre polynomial. In 
the derivation of Eq. (8), it has been assumed that the 
translational-rotational motions of a particle are decou- 
pled [34]. This decoupling is exactly valid for hydrody- 
namically interacting spheres at short times only. For 
non-spherical particles, it is an approximation even to 
linear order in t. 

At short times, Gr{t) = exp{— 6D^t} decays exponen- 
tially. At long times, however. Grit) decays in principle 
non-exponentially, with an average decay rate somewhat 
smaller than D^ [34]. While a genuine long-time rota- 
tional self-diffusion coefficient does not exist, one can de- 
fine instead a mean orientational self-diffusion coefficient, 
Dr, determined by the time dependence of Grit) for all 
times. A corresponding Dr of platelets is shown indeed 
in Fig. 2, as obtained in panel (d) of Fig. 1 using the 
KWW analysis. The calculation of Dr for colloidal hard 
spheres in [34] suggests that Dr is only slightly smaller 
than D^ , at least for smaller values of (p. The mean ro- 
tational diffusion coefficient D:^ depicted in Fig. 2, how- 
ever, decreases very strongly at larger (j), to an extent 
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Figure 4. Shear viscosity rj of gibbsite platelets in DMSO 
versus shear rate 7, for = 7.6%, 6.9%, 5.4%, 4.0%, 2.4%, 
and 1.45% (from top to bottom). 



comparable to that of Dt- Like for Dt, we attribute this 
strong decHne of D:^ at larger cj) to the strong anisotropic 
electro-hydrodynaniic coupling of the charged platelets. 



C. Static viscosity results 

It is interesting to note from Fig. 2 that the strong 
decline of the experimental Dt/D^ and Dr/D^ with in- 
creasing (j), is qualitatively similar to the concentration 
dependence of the normalized inverse low shear-rate lim- 
iting viscosity, rjo/rj. To obtain the low-shear rate vis- 
cosity data shown in Fig. 2, using an ArG2 rheometer 
from Ares, we have measured the viscosity as a func- 
tion of shear rate 7 for various concentrations. Results 
of these measurements are depicted in Fig. 4, not unex- 
pectedly showing that the gibbsite suspensions are shear- 
thinning, with the shear-thinning regime shifting to lower 
shear rates with increasing volume fraction. The low 
shear-rate limiting viscosity 77 is identified by the (mean) 
Newtonian plateau value at low shear rates, with results 
depicted by stars in Fig. 2. For technical reasons, no 
reliable viscosity data have been obtained for very low 
shear rates 7 < ls~^. Moreover, for technical reasons 
our viscosity data do not include the low-concentration 
regime from which the 0-independent intrinsic viscosity 
,[r]] = lim0^o('7 ~ '7o)/(?/o'?^')j could be determined, which 
depends on p only. For a sphere {p ~ 1) with stick hydro- 
dynamic boundary conditions, [77] = 2.5, and larger val- 
ues for [rj] are obtained with increasing asphericity. From 
the power-law representation of the simulation data in 
[10] for the intrinsic viscosity of thin cylinders, [77] w 6.8 
is obtained for p = 0.087. This value is similar to the one 
observed in [G-j] . Note that a spheroid of equal aspect ra- 
tio has a distinctly larger intrinsic viscosity of [77] « 9.1. 



IV. CONCLUSIONS 

Using (D)DLS and SLS, we have measured the long- 
time translational collective- and self-diffusion coeffi- 
cients Dc and D^, the mean rotational diffusion coef- 
ficient D^, and the static scattered intensity I{q) of 
charged gibbsite platelets in DMSO at low ionic strength. 
Our experiments cover the concentration range from very 
dilute systems up to the I/LC transition. 

Our usage of DMSO as a solvent with a dielectric con- 
stant close to that of gibbsite, has enabled us to deter- 
mine the translational and rotational diffusion properties 
without the necessity of invoking elaborate X-ray pho- 
ton correlation spectroscopy measurements. A fast re- 
laxation mode in the dynamic scattering data has been 
frequently reported in relation to the liquid-glass transi- 
tion, the glassy state and also for polymer coated clays 
[G, 8, G4]. In the present study, a fast mode has been 
found also for the isotropic phase, and it has been iden- 
tified as a collective diffusion mode. 

With increasing </>, the measured collective diffusion 
coefficient increases up to about twenty times the single- 
particle (orientationally averaged) translational diffusion 
coefficient. Different from the translational and rota- 
tional self-diffusion coefficients, which strongly decrease 
for (f) > 2%, Dc stays nearly constant for these larger 
concentrations even up to the I/LC transition concentra- 
tion. The strong decay of the self-diffusion coefficients 
is accompanied by a comparatively pronounced increase 
of the zero-shear-rate limiting static shear viscosity. We 
have provided arguments, both experimentally and the- 
oretically, that the coefficients Dt and D:^ obtained in 
our scattering modes analysis, should be identified, re- 
spectively, with the translational long-time and the mean 
rotational self-diffusion coefficients of gibbsite platelets. 

At low concentrations, I(q) is well reproduced in the 
simplifying translational-rotational decoupling method, 
where correlations between particle positions, sizes and 
orientations are neglected, and where the direct platelet 
interactions are approximated by a spherically symmetric 
electrosteric repulsion of DLVO type. The effect of hy- 
drodynamic interactions is accounted for in our effective 
sphere model. 

Except for very low concentrations, the accessible q- 
range in light scattering experiments is restricted to 
wavenumbers smaller than the value where the princi- 
pal structure peak in I(q) occurs. The effective particle 
charge has therefore been estimated, from the I/LC tran- 
sition concentration, as Z = 71, and was kept constant 
in our model calculations independent of concentration 
and salinity. For an unambiguous determination of Z, 
and to explore its dependence on cj) and ris, a broader 
g-rangc than accessible in SLS is required, which can be 
studied in future SAXS and XPCS measurements. In a 
more refined theoretical model for I{q), one can account 
for shape-dependent direct interaction contributions us- 
ing the PRISM model [G5-68]. However, regarding the 
dynamic quantities, it will be very difficult to include 
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shape-dependent hydrodynamic interactions, on avoid- 
ing numerically expensive multiparticle collision, fluctu- 
ating Lattice-Boltzmann or Stokesian dynamics simula- 
tions. We note here that while the simple effective sphere 
model clearly fails in terms of quantitative predictions, it 
allows for a correct assessment of qualitative features in 
the isotropic phase such as the ordering in magnitude of 
the rotational and translational self-diffusion coefficients, 
and the approximate plateau region of the collective dif- 
fusion coefficient at larger (p. 

Upon increasing the concentration, peculiar observa- 
tions have been made in the measured diffusion proper- 
ties of the gibbsite/DMSO system. Similar to suspen- 
sions of charge-stabilized colloidal spheres at low ionic 
strength, Dc grows with increasing concentration. Since 
the value of is basically the same at long and short 
time scales, and since the increase in Dc{<j}) is mainly 
due to the reducing osmotic compressibility S{q — > 0), 
these basic features of Dc{4>) have been reproduced in 
our effective sphere model, using the calculated S{q) as 
input to the otherwise parameter-free self-part corrected 
(57-scheme describing the short-time dynamics. 

The employed effective sphere model of gibbsite is 
less accurate regarding the measured long-time dynamic 
quantities and ZJ^f, which are more sensitive to the 
anisotropic direct and hydrodynamic interaction parts 



most influential on shorter length-scales. The strong de- 
cay of both self-diffusion coefficients to less than 1% of 
their respective infinite dilution values near (pi/i^c — 8% 
is not reproducible in an effective sphere model, where 
Dt ^ 0-1 X D^ even at the largest fluid-state concentra- 
tion. 

A concentration-dependence similar to the ones of 
and has been found for the inverse low shear-rate lim- 
iting viscosity {rj/rjt)) ^ . Examinations of possible gen- 
eralized Stokes-Einstein relations between viscosity and 
the various diffusion coefficients in concentrated platelet 
ffuids could be the topic of a future study based on the 
present work. 
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